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We investigate the impact of impurity scattering on superconductivity in an anisotropic multi¬ 
orbital model with spin-orbit coupling which describes the electron fluid at two-dimensional oxide 
interfaces. As the pairing mechanism is under debate, both conventional and unconventional super¬ 
conducting states are analyzed. We consider magnetic and nonmagnetic spin-dependent intra- and 
interorbital scattering and discuss possible microscopic realizations leading to these processes. It 
is found that, for magnetic disorder, the unconventional superconductor is protected against inter¬ 
band scattering and, thus, more robust than the conventional condensate. In case of nonmagnetic 
impurities, the conventional superconductor is protected as expected from the Anderson theorem 
and the critical scattering rate of the unconventional state is enhanced by a factor of four due to 
the spin-orbit coupling and anisotropic masses in oxide interfaces. 


I. INTRODUCTION 

To understand the properties of many of the recently 
discovered superconductors^’^, a multi-band description, 
which is often tied to the existence of multiple orbital 
degrees of freedom, is essential. Taking into account 
several Fermi surfaces does not only yield quantitative 
corrections, but also makes qualitatively different con¬ 
cepts possible such as pairing states that have no ana¬ 
logue in a single-band picture. One famous example is 
the 5“' state that has been proposed, e.g., in the iron- 
based superconductors^’'*^ and earlier in other systems®. 
While the pairing field is constant on a given Fermi sur¬ 
face, it changes sign between the electron and hole band. 

The identification of a pairing state and mechanism 
constitutes a challenging task, in particular, in a multi¬ 
band system. The behavior of superconductivity under 
the influence of different types of disorder can give impor¬ 
tant hints about the order parameter and pairing glue. In 
case of conventional s-wave superconductivity, the tran¬ 
sition temperature is only very weakly affected by the 
presence of nonmagnetic impurities except for very strong 
disorder where the mean-free path I is comparable to or 
smaller than the Fermi wavelength kp^. This is usually 
referred to as the Anderson theorem®^®. When magnetic 
impurities are considered, the transition temperature is 
suppressed® and vanishes when I becomes of the order 
of the superconducting coherence length ^ even when 
kpl ^ 1. Since various unconventional pairing states 
have been shown*^®^-*® to be strongly suppressed already 
by nonmagnetic impurities, the sensitivity of the transi¬ 
tion temperature to nonmagnetic disorder is commonly 
used as an indication of unconventional pairing. Un¬ 
conventional superconductivity is thus expected to occur 
only in very clean systems. 

However, several studies*^®’*^®^*^® of both magnetic and 
nonmagnetic disorder have revealed that, in some cases, 
unconventional superconductors can be much more ro¬ 
bust than naively expected or even enjoy an analogue 
of the Anderson theorem. Refs. 12 and 16-18 indicate 


that it is essential to take into account all characteris¬ 
tic details of the system, such as multiple orbitals and 
spin-orbit coupling, in order to judge reliably the robust¬ 
ness of unconventional superconductors in the presence 
of disorder. From a theoretical point of view, it then be¬ 
comes more difficult to And an approach that captures 
the essential features but leads to results that are uni¬ 
versal in the sense that they do not depend on irrelevant 
microscopic details. 

One example of a recently discovered superconducting® 
system where the multi-orbital character*^®’®®, spin-orbit 
coupling®*^’®® as well as strongly anisotropic effective 
masses®® seem to be crucial for understanding its prop¬ 
erties is provided by the two-dimensional (2D) electron 
liquid that forms®'* at the interface between the insulators 
LaA102 (LAO) and SrTi02 (STO). In Ref. 25, two possi¬ 
ble pairing states have been proposed for this system. If 
superconductivity arises as a consequence of conventional 
electron-phonon coupling, the pairing held will have the 
same sign on both spin-orbit split Fermi surfaces. In 
case of an unconventional, i.e., purely electronic, pairing 
mechanism, the order parameter has opposite signs on 
the two Fermi surfaces. Despite the similarities to the 
s“* state known from the iron-based superconductors, 
there is one main crucial difference. The Fermi surfaces 
in LAO/STO interfaces are singly degenerate due to the 
combination of broken inversion symmetry and spin-orbit 
coupling which has important consequences for supercon¬ 
ductivity (see also, e.g.. Refs. 26 and 27). For example, 
the unconventional pairing state proposed®® for oxide in¬ 
terfaces has been shown to be a time-reversal invariant 
topological superconductor which is intimately related to 
the Fermi surfaces being split by spin-orbit coupling. 

Consequently, it constitutes an important open ques¬ 
tion whether such an exciting new state of matter is in¬ 
deed realized in oxide interfaces. An analysis of the im¬ 
pact of disorder on the superconducting states proposed 
in Ref. 25 is thus required. It is particularly impor¬ 
tant to perform a quantitative analysis beyond a simple 
one-orbital toy model calculation since experimental es- 
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timates of the coherence length^ ^ ~ 70 —105 nm and the 
mean-free path^®’^® / ~ 25 — 95 nm are of the same order. 

In this paper, we consider impurity scattering in super¬ 
conducting LAO/STO interfaces using a two-band model 
that captures the multi-orbital character, the spin-orbit 
splitting as well as the anisotropy of the masses. The 
resulting highly non-circular form of the Fermi surfaces 
and anisotropic textures of the wavefunctions are con¬ 
veniently taken into account by introducing patches in 
the nested regions of the Fermi surfaces. Due to the 
strong anisotropy, the wavefunctions are approximately 
constant within each patch yielding results that only de¬ 
pend on a small set of parameters. We find that the spin 
and orbital polarization of the states at the Fermi sur¬ 
faces enhances the stability of both the conventional and 
unconventional superconducting phase against impurity 
scattering. Furthermore, contrary to common wisdom, 
the unconventional superconductor is found to be, by a 
factor of two, more robust against magnetic scattering 
than the conventional condensate. 

The remainder of the paper is organized as follows. In 
Sec. II, we introduce the clean model of superconduc¬ 
tivity. The impurity configurations consistent with the 
symmetries of the system and possible microscopic real¬ 
izations thereof are discussed in Sec. III. In Sec. IV, we 
present the approximations used to calculate the transi¬ 
tion temperature in the disordered system. Finally, the 
results can be found in Sec. V. 


II. SUPERCONDUCTIVITY IN THE CLEAN 
SYSTEM 

To describe superconductivity in oxide heterostruc¬ 
tures, we use an effective Hamiltonian of the Ti id^z 
and ‘idyz orbitals. To justify why these are the essen¬ 
tial degrees of freedom, let us recall that the orbitals 
relevant for the electronic low-energy physics are the 
Ti Zd states of the t 2 g manifold, i.e. the Zd^y, Zd^z 
and Zdyz orbitals. The states of the interface elec¬ 
tron system derived from the Zdxy orbital mainly lo¬ 
calize and possibly order magnetically®®^®®. This ex¬ 
plains the discrepancy between the carrier concentration 
in the two-dimensional electron liquid expected from the 
polar catastrophe mechanism and that seen in trans¬ 
port measurements. The electronic states that lead to 
superconductivity are therefore most likely due to the 
Zdxz and Zdyz orbitals which is supported by recent 
experiments®®’®®. 

Let us first focus on the non-interacting part of the 
effective Hamiltonian of these states, 

(1) 

k 

where il)a,k and f, describe the annihilation and cre¬ 
ation of quasi-particles of crystal momentum k in state a, 
with a labeling the four distinct combinations of the 
two orbitals {Zdxz, Sdyz) and the spin orientations. The 


Hamiltonian is characterized by two ingredients, h{k) = 
^m(^) T hso{^)- 

Firstly, the different overlap of the two orbitals along 
the X- and y-direction leads to anisotropic masses de¬ 
scribed by 
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where the upper (lower) component refers to the Zdxz 
{Zdyz) orbital and ctq is the unit matrix in spin space. Ex¬ 
perimentally, large anisotropy ratios m\ilm\ ~ 15... 30 
have been reported®®. The orbital mixing term 5 de¬ 
scribes interorbital second-nearest neighbor hopping. 

Secondly, the properties of the electron liquid are 
strongly affected by spin-orbit coupling®®’®®. As shown 
in Ref. 25, this leads to the contribution 


hso{k) = -/3t 2(T3 -I- QfiTo (fciCT 2 - k 2 (Ji) 

a2Ti (fciCTi - k2Cr2) + asTs {ki(T2 + fecTi) 


( 3 ) 


to the Hamiltonian. Here ai and (3 are real constants and 
CTi {xi) denote Pauli matrices in spin (orbital) space. We 
see that besides the usual Rashba coupling (oc ax), the 
multi-orbital character of the model also allows for the 
atomic spin-orbit term proportional to /3 and two Dres- 
selhaus couplings ( 02 , 03 ). Eq. (3) is the most general 
spin-orbit coupling Hamiltonian (up to linear order in k) 
consistent with time-reversal and the C^y point symme¬ 
try of the system. Starting from a three-orbital model 
that also includes the Zdxy band and integrating out®® 
the latter reproduces Eq. (3) with the additional con¬ 
straint ax = —02 = — 03 - 

Note that we do not include any coupling to the mag¬ 
netic moments that possibly develop in the Zdxy bands. 
The reason is that no significant correlation between the 
inhomogeneous magnetization of the Zdxy orbitals and 
the superconducting response has been observed®^ and, 
hence, the effective exchange coupling between the local¬ 
ized magnetic states and the superconducting electrons is 
expected to be negligible for describing superconductiv¬ 
ity. This can be understood as a consequence of the spa¬ 
tial separation of the localized magnetic and the itinerant 
superconducting bands and the small orbital admixture 
of Zdxy in the states that host superconductivity. 

The spin-orbit coupling (3) breaks inversion symmetry 
and removes the spin-degeneracy of the Fermi surfaces. 
For chemical potentials close (deviation < /3) to the bot¬ 
tom of the spectrum of h, only two out of the four singly 
degenerate bands of Eq. (1) are relevant. As carrier con¬ 
centration dependent studies®®’®® show that the sweet 
spot of superconductivity is associated with the Fermi 
level entering the xz- and j/z-bands, we will focus on the 
two bands of h{k) with lower energy. The corresponding 
form, spin and orbital polarization of the Fermi surfaces 
are shown in Fig. 1. In the 16 patches, highlighted in 
green, the Fermi surfaces are nearly straight lines which 


3 


{xz) 



(xz) 



Figure 1. Outer (a) and inner (b) Fermi surface (black lines) 
with the corresponding spin (red arrows) and orbital (color) 
polarizations within the two-orbital model (1). The patches 
with approximately constant wavefunctions are indicated in 
green. Here we have used /3 = 20meV, ai = lOmeVA, 
mh/mi = 30, mi = 0.7me (me is the free electron mass) 
and 5 — 0.28m^^ deduced from Refs. 19, 21, 23, and 35. 


are, in groups of four, approximately parallel (nesting). 
Within each patch, the wavefunctions are almost con¬ 
stant. The large mass anisotropy in leads to strong 
orbital polarization: The patches in the lobes along ki 
(^ 2 ) are mainly j/z-polarized (ccz-polarized). Similarly, 
the spin-orbit coupling hgo largely aligns the spin parallel 
or antiparallel to one of the in-plane coordinate axes. The 
strong spin and orbital polarization of the patch wave- 
functions is of crucial relevance for impurity scattering 
as it can lead to the suppression of certain processes. 
This will be discussed in detail in Sec. V. 

In Ref. 25, it has been shown that the system can have 
two distinct superconducting instabilities: In case of a 
conventional, electron-phonon induced, pairing mecha¬ 
nism, the pairing field has the same sign on both Fermi 


surfaces, which will be referred to as s’*"’*' in the follow¬ 
ing. If the superconductor is due to the exchange of 
particle-hole fluctuations, the sign of the order parame¬ 
ter will change between the Fermi surfaces. Correspond¬ 
ingly, this unconventional superconducting state will be 
denoted by . 

At low energies, the relevant interaction channel that 
distinguishes between these two superconductors is pair 
hopping, i.e., the annihilation of a Kramers pair of elec¬ 
trons from one band and the creation of a Kramers pair 
in the other. Depending on the sign of the associated 
coupling constant, either the or s"' superconductor 
will be favored. As the sign conventions for this interac¬ 
tion term sensitively depends on the phase convention for 
the eigenstates of h, the mathematical formulation of the 
pair hopping interaction will be postponed to Sec. IV B. 


III. DISORDER CONFIGURATIONS 


Here we discuss the possible impurity configurations 
that can be present in oxide interfaces. The most gen¬ 
eral (quadratic) Hamiltonian of a given disorder realiza¬ 
tion W reads 


Hdis = f (4) 

J x,x' 

where, for convenience, we have switched to a real- 
space representation with ipaix) and denoting the 

Fourier-conjugates of V'a.fc and '0^ respectively. 

In what follows, we not only allow for disorder config¬ 
urations that are consistent with the point group of the 
system but also take into account impurities that locally 
break any subset of the point symmetries. We only re¬ 
quire these symmetries to be preserved on average. For 
example, oxygen vacancies that locally break the fourfold 
rotation symmetry occur with equal probability on bonds 
along the two perpendicular directions parallel to the in¬ 
terface. This effect is most easily discussed on the level of 
correlation functions. Focusing on Gaussian-distributed 
disorder, the entire information about the statistics is 
contained in the correlator 

Baja' 3:2) 

(5) 

:= (WAia; ixi,x[)Wa^a'^ ix2, x' 2 ))^.^. 

Here (.. .)^;g represents the disorder average only tak¬ 
ing into account Hermitian disorder configurations, W = 
FF^. The latter constraint makes the averaging procedure 
equivalent to a Gaussian integration over a real-valued 
field. 

Within the replica approach^®, averaging over W leads 
to an effective impurity-induced electron-electron inter- 












4 


action that reads in the action description as 


the approximation (7), to the constraint 


^ r.r' = l'^ J ^1,^1,^2, 31 2 

X (a;i, (x2, t')C' K. 't') 

^ q, 2 ck' ( 3 ?!, 3 ?]^, ^^ 2 , ^^ 2 ). 


( 6 ) 


Here, ip'^, ip^ are the Grassmann analognes of ipa, V’a 
with the additional replica index r = 1,..., i?, where R 
denotes the nnmber of replicas. At the end of the calcu¬ 
lation of physical quantities, the limit i? —?► 0 has to be 
taken. We see that the correlator defined in Eq. (5) plays 
the role of the bare disorder vertex function. The same 
holds in the diagrammatic disorder-averaging technique 
of Ref. 37. 

For simplicity, we will restrict our analysis to spatially 
local disorder, i.e. Xj = 3 ;' in Eq. (6), and ^-correlated 
statistics which, in addition, forces Xi = x^- Further¬ 
more, let us assume the disorder distribution to be ho¬ 
mogeneous. Taken together, the impurity vertex becomes 
spatially local and translation invariant. 


,a2CK2 (^1’ ^1: ^2, ^ 2 ) 

= d{xi - x[)S{x 2 - x' 2 )d{xi - 3;2)r;^„, . 


(7) 


A. Symmetry properties of the impurity vertex 

To derive and classify the disorder configurations rel¬ 
evant for oxide interfaces, we have to analyze how the 
vertex T^^^, transforms under the point group op¬ 
erations and time-reversal. 

To begin with the former, consider an arbitrary oper¬ 
ation g of the point group. The field operators transform 
according to 

ipa{.x) 4 (nUg)) ijjp(TZyig)x), (8a) 

ipiix) 4 tpl{TZy{g)x){TZ^{g))^^, (8b) 

where Tilp{g)JTZji,{g) = TZs{g)J with TZvig) and Ti-sig) 
denoting the representation of g on vectors and pseu¬ 
dovectors, respectively. Here J represents the total an¬ 
gular momentum operator. From Eq. (4) readily follows 
that 


r: 


QllQj^ , 0 : 20:2 






( 10 ) 


for all operations g. 

To distinguish between nonmagnetic and magnetic dis¬ 
order, we have to take into account the properties under 
time reversal. Let us denote the representation of time 
reversal on wavefunctions by 0 = TA with unitary R 
and K. representing complex conjugation. Within the 
basis used in Eq. (3), it holds T = iToa 2 - According 
to Eq. (4), the disorder Hamiltonian transforms under 
time-reversal as 

W(3:,3;') 4 TW^{x',x)T-^. (11) 


Here we have already taken advantage of the Hermiticity 
of W which allows for a representation that does not ex¬ 
plicitly involve complex conjugation. This is convenient 
as the restriction on time-reversal symmetric (TRS), 
CT = 1, or time-reversal antisymmetric (TRA), a = —1, 
disorder realizations can then be reformulated in terms 
of the constraints 


r' 

,OL20l'2 


^'Taicti^a[ai,01202 )ajaj 

C’’42 02^0 id'j,5^52 ) a’a '2 


( 12 ) 


on the disorder vertex. Note that it is sufficient to impose 
the condition in the first line as the second line automat¬ 
ically follows from the property 


^aiaj , 020^2 ^020*2^010'^^ ( 73 ) 

which is a direct consequence of the definition (5) . 

Because of the very different behavior of superconduc¬ 
tivity in the two cases, we will discuss TRS and TRA 
disorder separately, which will be referred to as “non¬ 
magnetic” and “magnetic” in the following. 


B. Nonmagnetic disorder 

To write down the most general disorder vertex 
^'aia[, 02 a '2 Satisfying Eqs. (10) and (12) with cr = 1, 
we perform an expansion in terms of Hermitian matri¬ 
ces {w^}, i.e., let 


W(3;,3;') 4 n^{g)W {np\g)x,Tl-\g)x')Tllig), 

(9) 

which immediately yields the transformation behavior of 
the correlator (5). As mentioned, while a generic given 
disorder realization W will break all point symmetries of 
the system, on average the full point symmetry of the 
normal state has to be restored. In other words, the cor¬ 
relator (5) must transform trivially which leads, within 


^aia[,a2a’2 ~ {w^i)^ 202 - (74) 

Note that as a consequence of 

= W in Eq. (5). Together with Eq. (13) this forces 
C to be real and symmetric, 

C = C* = C^. 


(75) 
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TABLE I. The symmetry properties of the basis functions 
constructed from the Pauli matrices in spin [oi) and orbital 
(Ti) space. Here Ai, A 2 , B\, B 2 and E denote the irreducible 
representations of C' 4 „, whereas TRS (TRA) indicates that 
the matrix is even (odd) under time-reversal. 



To 

Tl, 02 

03 

To 

Ai/TRS 

E/TRA 

Aa/TRA 

n 

R 2 /TRS 

E/TRA 

Hi/TRA 

r2 

Aa/TRA 

E/TRS 

Ai/TRS 

T3 

Hi/TRS 

E/TRA 

Ha/TRA 


Let us represent the basis matrices as the 16 

possible tensor products of the Pauli matrices Ui and Ti, 
i = 0,1,2, 3, referring to the spin and orbital degrees 
of freedom as in Eq. (3). This constitutes a convenient 
choice since these basis matrices both transform under 
the irreducible representations of the point group 
and are either symmetric or antisymmetric under time- 
reversal as summarized in Table I. Recalling Eq. (15), 
one can then readily construct the most general vertex 
in case of nonmagnetic disorder, 


r' , 

ot\(y.-^,cx.2C>L2 


- lAi 


+ 




(tzCTs) 

c 

+ 7Si ('r3cro)ai„' (racro) 


+ Tl 


0 : 20:2 

(T20-l)„ja( ('^20 -i)„2„^ 


+ 7 Ai ('^ 20'3)a2Q,^ 

^ (:^0O'0)a2a( 

+ IB2 (Ticro)„ic; 0202 

+ (r2cr2)„i„- (r2Cr2) 

0202 J 

(16) 


characterized by six independent real coupling constants. 
To gain physical insight into Eq. (16), let us first discuss 
the meaning of the different terms and then analyze how 
they can emerge microscopically. 

The contributions proportional to 75 ^ and 
correspond to spin-trivial disorder with orbital wavefunc- 
tions of different point symmetry: The different terms 
transform as Ai, Bi and B 2 under which corresponds 
to S-, dx2-y2-, and d^y-wave. Noting that T2 is the projec¬ 
tion of the z-component of the orbital angular momentum 
operator on the subspace spanned by the 3dxz and 3dyz 
orbitals, the terms with 7 ;^^ and je represent (atomic) 
spin-orbit disorder with out-of-plane/longitudinal and in¬ 
plane/transversal orientation of the spin, respectively. 
The prefactors of the two terms in the last line of Eq. (16) 
must be identical in order to ensure C' 4 -rotation symme¬ 
try of the disorder distribution. Finally, the contribution 
proportional to 7 )^^ describes the fact that the C' 4 „ point 
group allows for an admixture of s-wave disorder and 
longitudinal spin-orbit disorder on the same site. 

To understand how the different terms can be realized 
microscopically let us investigate the schematic illustra¬ 
tion of the crystal structure in Fig. 2. We begin with 




Figure 2. Schematic of four unit cells of the crystal structure 
of STO. The part shown is assumed to be close to the in¬ 
terface (parallel to the a:?/-plane) in the spatial region where 
the conducting electron system resides. The numbers refer to 
the different positions of Ti (1), O (2, 3, 3') and Sr (4, 4') 
which are referred to in the main text for the discussion of 
the different terms in the impurity vertices. 


oxygen vacancies, which are frequently discussed in the 
context of LAO/STO interfaces (see, e.g., Refs. 38 and 
39). Assume that one oxygen at position 2 , right be¬ 
low the central Ti at position 1, is vacant. In this case, 
the defect preserves the full C 4 y point symmetry of the 
crystal such that only the first 3 terms in Eq. (16) can 
be present. The explicit values of 7 ) 4 ^, 7 ) 4 ^ and 7 ;^^ de¬ 
pend on microscopic details of the system, however, as 
will be seen in Sec. V, only the sum of these scattering 
rates enters the effective pair-breaking strength. Here we 
mention that all coupling constants in the impurity ver¬ 
tex are proportional to the concentration of the defects 
considered. If instead an oxygen adjacent to the central 
Ti, e.g., at position 3, is vacant, all point symmetries 
of the lattice are locally broken except for the mirror 
symmetry with respect to the j/z-plane. This symmetry 
constraint rules out the contribution proportional to 73 ^ 
and the first term in the last line of Eq. (16) as both tiCTo 
and T 2 (Ti are odd under {x,y) —>■ {—x,y). Correspond¬ 
ingly, if we consider oxygen vacancies at position 3', the 
second term in the last line of Eq. (16) and again the 
contribution with 75 ^ will be forbidden. Recalling that 
the prefactors of the two terms in the last line of Eq. (16) 
have to be equal by C 4 rotation symmetry, we see that, 
as expected, vacancies at position 3 and 3' have to be 
included with equal probability to restore the full point 
symmetry after averaging. In summary, single oxygen 
vacancies allow for all contributions in Eq. (16) except 
for the 7 _b 2 “ term. 
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To see how the latter can occur, let us now assume that 
a Sr atom is vacant, e.g., at position 4 in Fig. 2. Similarly 
as above, all point symmetries are broken save for the 
mirror reflection at the (110) plane. This now renders 
the 7 B 2 “term possible which is even under the reflection 
at (110). Assuming that, with the same probability, a 
Sr atom is missing at position 4', one restores rotation 
symmetry. In this case, all terms in Eq. (16) are allowed 
except for the one proportional to 75 ^. Finally, we note 
that the same is true for a local oxygen double vacancy 
where oxygen atoms are missing at position 3 and 3' at 
the same time. We have thus provided examples of how 
all the terms in the most general non-magnetic impurity 
vertex (16) can be realized. 


C. Magnetic disorder 


In case of magnetic disorder, i.e. a = —1 in Eq. (12), 
we proceed in the same way as in Sec. IIIB. From Ta¬ 
ble I, one finds that there are three independent terms in 
the impurity vertex that correspond to in-plane magnetic 
moments. 


p/ II 

CXiCx'-^ ,Q:2Ck2 


ik 

2 

2 

zi 

2 


('^l<^l)aia; i'^i^i)a2a'2 + a2a'2 




(Tt^l)cK2a2 (^5^2)^^^^ (Ti^2)Q,2a2 

(17a) 


and five terms. 


F' 


CXiCX-^ , 02(^2 


= 7^2 iTo<y 3 )a,a[ ('^00-3)0,20^ + 7 Bi (ncra)^^^/^ (TlCra)^^^^ 

+ IB 2 iT3Cr3)a2a'2 + 7^2 




(T2Cro)„iai No-3)„2ai + ^ 

(17b) 


where the magnetic moment is oriented perpendicular to 
the plane. As before, all coupling constants must be real 
due to Eq. (15). 

Physically, the three terms 7 ^, 7 ^ and 7 )^^ in F'^' cor¬ 
respond to impurities with in-plane spin-magnetization 
and distinct orbital symmetry (s-, dxy- and dx 2 -y 2 -'wa,ve). 
Note that no orbital-magnetic disorder with in-plane ori¬ 
entation can occur in the two-orbital model as the pro¬ 
jection of the X- and y-component of the orbital angular 
momentum operator onto the 3dxz and Sdyz subspace 
vanishes identically. The analogous terms in r'"*“ with 
the same orbital symmetry but spin polarization along 
the z-direction are proportional to JBi and 7 ^ 2 ) 
respectively. Finally, 7 )^^ corresponds to purely orbital 


magnetism along the z-direction and 7 ™ shows that an 
arbitrary admixture of spin and orbital magnetism along 
the z-direction is allowed by symmetry. 

Similarly to the nonmagnetic case, let us discuss some 
examples of how the terms in Eq. (17) can emerge mi¬ 
croscopically. One can imagine that locally one Ti atom 
orders magnetically due to the proximity of the system 
to a competing magnetic instability^®. Numerical calcu¬ 
lations indicate"^^’'^^ that magnetic moments on Ti atoms 
can be induced by the presence of oxygen vacancies. Al¬ 
ternatively, we may think of Ti being replaced by a differ¬ 
ent magnetic atom. In all cases, the form of the resulting 
impurity potential crucially depends on the orbitals that 
host the magnetic moment and whether it is orbital or 
spin magnetism. Thus, in general, all terms in Eq. (17) 
are feasible. To be more specific, recalling that the ob¬ 
served magnetism at the interface is mainly due to the 
3dxy band®°^®®, one may expect that, most likely, the spin 
degree of freedom of the Ti 3dxy orbital orders locally. In 
that case, the impurity vertex would be described by the 
terms proportional 7 ^ and 7^1 as the associated orbital 
structure, ri, transforms as xy under C^y. 

If, instead, a magnetic moment develops at a site that 
does not lie at a high symmetry point with respect to 
the Ti atoms, the symmetry constraints will be less re¬ 
strictive. Assuming, e.g., a spin-magnetic moment along 
the z-axis with s-wave orbital symmetry (i.e. tocts) at 
the oxygen site 3 (3') in Fig. 2, the disorder potential is 
only restricted to be odd both under time-reversal and 
mirror reflection symmetry with respect to the j/z-plane 
(a;z-plane). The reduction of symmetry constraints sim¬ 
ply follows from the fact that, irrespective of the sym¬ 
metries of the local impurity, the presence of a pertur¬ 
bation at 3 (3') already breaks all symmetries except for 
time-reversal and one mirror reflection. Under these two 
remaining symmetry operations, the resulting impurity 
potential in the effective model for the conducting Ti or¬ 
bitals must have the same behavior as toCTs (in both cases 
being odd). From these criteria one finds that only the 
term proportional to 75 ^ can be excluded in Eq. (17). 


IV. CRITICAL TEMPERATURE OF THE 
DISORDERED SUPERCONDUCTOR 

In this section, we provide details about how the crit¬ 
ical temperature of the superconductors is calculated in 
the presence of disorder. Readers that are mainly inter¬ 
ested in the results of this paper and the physical inter¬ 
pretation thereof can skip this part and directly proceed 
with Sec. V. 


A. Patch approximation 

For our calculation, it will be convenient to work in 
the eigenbasis of the normal state Hamiltonian h{k) in 
Eq. (1). For this purpose, we introduce new operators 
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fa,k and fl ,^ via 

= (0/3,fc)„ //3,fc, ^i^k = fl,k > (18) 

where (j)a,k are the four eigenstates of h{k). By design, 
the normal state Hamiltonian (1) becomes diagonal, 

= (19) 

k 

In the new basis, the momentum space representation of 
the impurity vertex in Eq. (7) reads 

{4*011,ki) ai {4*012,k-^ 02 1^3i5j,52S2 

( 20 ) 

which depends on all four external momenta. To sim¬ 
plify the analysis, we take advantage of the fact that, 
within each of the patches of the Fermi surfaces in Fig. 1, 
the wavefunctions (po.k are approximately constant. To¬ 
gether with the fact that this set of patches covers the 
main part of both Fermi surfaces, it is justified to replace 

4^o,k ^ fo,k ^ /(Aj,r/),p? 

( 21 ) 

a,fe XJ,ri,p 

We only keep the fermions from the patches as dynam¬ 
ical degrees of freedom characterized by a constant spin 
and orbital polarization in a given path. This is what 
we refer to as the “patch approximation”. In Eq. (21) 
and in the following, A = 1, 2 denotes the two Fermi sur¬ 
faces and j = 1, 2, 3,4, 77 = -I- (77 = —) represent the four 
distinct patches on each Fermi surface in the half-space 
with ki + k 2 >0 (fci -I- ^2 < 0) as illustrated in Fig. 3. 
The additional index 77 will be convenient for describing 
superconductivity as, for fixed A and j, the states with 
77 = -|- and t] = — are related by time reversal. We use 
p to denote the deviation of the crystal momentum from 
the centers of the different patches. We always assume 
that these momenta are cut off in a way that does not al¬ 
low for patches to overlap. Note that applying the patch 
approximation automatically takes advantage of the fact 
that, as discussed in Sec. II, two out of the four bands 
of the normal state Hamiltonian (1) can be neglected for 
energetic reasons. 

Using the compact multi-index notation t = (A,j, 77 ), 
the disorder induced interaction ( 6 ) finally becomes 
within the patch approximation 



5'r[V',V'] = -o 


E 


r 


PA 


Pl,Pl,P2 7P2 




P2P2’ 

( 22 ) 



Figure 3. Illustration of the different patches (green shaded 
domains) on the Fermi surfaces and the notation we use to 
label them. The patches are indicated by a composite index 
T = (A, j, 77 ) with A, j, and 77 referring to the two Fermi sur¬ 
faces, the four patches with fixed sign of fci + k 2 , and to the 
two half-spaces, respectively. 


where, using the expansion (14), 



In the second line of Eq. (23), we have introduced the 
diagrammatic symbol for the impurity vertex that will 
be used in Sec. IV C. Due to the translation invariance of 
the disorder distribution (see Eq. (7)), the impurity in¬ 
duced interaction (22) must conserve momentum, which 

is accounted for by G |0,1|. Here we will not 

specify the explicit form of D. Let us only note that it 
crucially depends on details of the band structure. To see 
this, consider, e.g., the scattering process of an electron 
from (1, 2, -b) to (2, 2, +) with an electron that starts and 
does not leave the patch (1,4,—) in Fig. 3. Obviously, 
the available phase space of this process depends both 
on the nesting vector between (1,2,4-) and (2,2,4-) as 
well as on the geometry of the patches. Fortunately, as 
we will explain in Sec. IV C, the transition temperature 
of the superconducting states will not depend on these 
details of the band structure. 

We emphasize that the vertex T^^ contains the full 
information about disorder scattering within the patch 
approximation, combining knowledge about the micro¬ 
scopic structure (see Sec. HI) of the impurities involved 




and the form of the eigenstates (see Fig. 1) between which 
electrons are scattered. Motivated by the strong orbital 
and spin polarization of the patch wavefunctions, we ap¬ 
proximate the states to be perfectly polarized. More 
precisely, within each patch, we take the spin to be fully 
aligned parallel to one of the coordinate axes and the or¬ 
bital weight to be either fully xz or yz corresponding to 
the approximate polarization of the wavefunctions illus¬ 
trated in Fig. 1. This simplification allows us to obtain 
analytic expressions for the scattering vertex F^^. In 
Sec. V, we will also comment on the impact of deviations 
from perfect polarization. 

B. Phase convention for the eigenstates 

To have a well-defined representation of the supercon¬ 
ducting order parameter in the eigenbasis of the normal- 
state Hamiltonian, it is essential to constrain the phases 
of the eigenstates (pa^k in Eq. (18). To see this, consider 
the microscopic order parameter Aapik) defined by the 
mean-field Hamiltonian 

^ E {'^lk^Mk)^l,-k + H.c.) , (24) 

k 


parameter transforms according to (again switching to 
the patch notation) 

VV' (28) 

where denotes the patch obtained when apply¬ 

ing the point group operation g on the subset of momen¬ 
tum space belonging to patch (A,j, 77 ). 

General symmetry and energetic arguments^^ indi¬ 
cate that the superconducting order parameter of ox¬ 
ide interfaces transforms according to a one-dimensional 
irreducible representation of the point group 
Experiments'^^ indicate that the superconductor is fully 
gapped. Together with the result that time-reversal sym¬ 
metry is not broken^^, the order parameter must trans¬ 
form under the trivial representation of C^y (s-wave) and, 
hence, 

^(Aj,? 7 ) = Aa. (29) 

The microscopic calculation of Ref. 25 yields Ai = A 2 
and Ai = — A 2 corresponding to the and s“' super¬ 
conductors as possible candidates. As already discussed 
in Sec. H, the crucial interaction leading to these two 
instabilities is pair-hopping, which is given by 


which reads in the eigenbasis (18) as 

A^,k ■■= (25) 

Here it has already been taken into account that we fo¬ 
cus on order parameters that have no Fermi-surface off- 
diagonal elements^^. Clearly, Eq. (25) is not invariant 
under a general fc-space local C/(l)-gauge transformation 
<t>a,k —>■ e.^^'^4>a,k of the wavefunctions. 

To remove this ambiguity, we take advantage of the 
fact that the normal state Hamiltonian h{k) is even under 
time reversal. Since the Fermi surfaces are singly degen¬ 
erate, one can always adjust the phases of the eigenstates 
to satisfy 


Hint-V E 

ij'=i p,p' \ 

^ /(2j',-),-p'/(2,i',-|-),p' +H.C. 


(30) 


Within the phase convention (26), taking V < Q {V > 
0 ) leads to the conventional (unconventional ) 
superconductor. 


C. Ginzburg-Landau expansion in the presence of 
disorder 


*^^(Aj', —r;)) ( 26 ) 

where, for convenience, we have switched to the patch 
notation. The phase factor —i is actually arbitrary but 
has been chosen to reproduce the conventions of Ref. 25. 
Note that the sign factor 77 is necessary due to the prop¬ 
erty 0^ = —1 of the time-reversal operator. Imposing 
Eq. (26) fixes the relative phase between (j)a.,k and 
thus, rendering Eq. (25) gauge invariant. 

To classify the possible superconducting phases, it is 
essential to know how the order parameter transforms 
under the elements of the point group of the system. Ap¬ 
plying Eq. ( 8 ) in Eq. (24) one readily finds that 

A(fc) 4 n^{g)A{TZp\g)k)TZlig). (27) 

Taking into account the phase convention (26) and using 
one more time that the Fermi surfaces are singly degen¬ 
erate, it follows that the representation (25) of the order 


To determine the transition temperature of the su¬ 
perconductor, we first decouple the pair-hopping in¬ 
teraction (30) by a mean-field approximation in the 
Cooper channel assuming s-wave symmetry as defined 
in Eq. (29). This yields 


^r = ^EE 


X,j,V P 


+ R.c.)-y{A*A,+AlA,). 


(31) 


Let us denote the free energy for a given disorder real¬ 
ization W in Eq. (4) by J^[W,Aa]. Assuming that the 
system is self-averaging, the free energy of the supercon¬ 
ductor is given by the disorder average AA])(jis that 

only depends on the disorder correlator (5) and not on 
details of the microscopic W. As we are only interested 
in the transition temperature of superconductivity, it is 
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Figure 4. As shown in (a), the coefficients dxx' of the 
Ginzburg-Landau expansion (32) can be represented using the 
full Greens function (double fine) and the dressed vertex (gray 
triangle). Part (b) and (c) illustrate the diagrams contribut¬ 
ing to the Greens function and vertex that are leading in the 
limit of weak disorder kpl ';$> 1. 


sufficient to expand the free energy up to second order in 
the order parameters. One finds 

A,]),;, ~ Aldxy{T)Ax, - ^ {A*A, + A^A^), 

(32) 

as Aa —>■ 0 , where the coefficients are given by 

^/(A,i.-r)),-p('r)/(Aj,7)),p('r) 

3,V,P 

X y]] V , 3 ' j' -ri')-p'{T )j J 

j'.v'-.P' /o/dis 

(33) 

in the path integral representation. Here (.. .)g denotes 
the expectation value with respect to the action associ¬ 
ated with the normal state Hamiltonian (19). The tran¬ 
sition temperature and the order parameter then simply 
follow from diagonalizing the quadratic form in Eq. (32). 
These two quantities are determined by the eigenvalue 
that first changes sign and the associated eigenvector, 
respectively. 

The disorder averaging in Eq. (33) leads to the effec¬ 
tive interaction (22) with the impurity line as defined in 
Eq. (23). The exact expression for the coefficient d\y{T) 
is represented diagrammatically in Fig. 4(a) and requires 
knowing the full self energy and vertex correction. In 
this paper, we will assume that the disorder scattering is 
sufficiently weak in the sense that the mean-free path I 
is much larger than the inverse Fermi momentum kp^. 
This assumption is justified by the experimental obser¬ 
vation that the typical normal state sheet resistance is 
by more than one order of magnitude smaller than h/e^, 
where e and h denote the electron charge and Planck’s 
constant.Being suppressed by a factor we 

neglect all diagrams with crossed impurity lines in the 
self-energy as well as in the vertex correction. Conse¬ 




quently, only the “rainbow diagrams” and “Cooperon lad¬ 
der”, which are shown in Fig. 4(b) and (c), have to be 
taken into account. 

From the diagrammatic representation, one observes 
that the impurity vertex T^^, , does not enter as a 

full tensor but only in form of the index combinations 


and 


jtt ^ rr' ,t't 


nPA 


- r'r,T' 


T ,P -►1 

_ 1 

~ 1 

- ► T, p 

r, p —► * 

— ► t ' ,p' 

r,p —► 

1 —► r',p' 

1 


i" 

1 


(34) 


(35) 


in the self-energy and vertex correction, respectively. 
Here the patch momenta p and p' have been added 
for later convenience and tk = (A, j, —rf) denotes the 
Kramers partner of patch r = (X^j,r]). In Eq. (35), we 
have taken into account that the superconducting ver¬ 
tex only couples Kramers partners, {r, p;tk,— p}, char¬ 
acterized by zero total momentum. As the patches are 
assumed to be disjoint in momentum space, momentum 
conservation allows for Cooper pairs only being scattered 
into pairs of states of the form {t' ,p'; —p'}. 

Note that this crucially reduces the complexity of im¬ 
purity scattering as the tensor T^'^ with 16^ possible in¬ 
dex combinations is replaced by two scattering matrices 
S and V with only 16^ entries each. More importantly, 
we directly see that the momenta p and p' can freely 
run over all values within the associated patches since 
momentum conservation is fulfilled in the diagrams of 
Eqs. (34) and (35) for arbitrary and independent values 
of p and p'. We emphasize that this crucial simplifica¬ 
tion is a consequence of the superconducting state be¬ 
ing characterized by Cooper pairs with vanishing center 
of mass momentum and the absence of diagrams with 
crossed impurity lines in the diagrammatic expansion of 
Fig. 4. Had we included diagrams with crossed lines, 
the impurity vertex would have entered in a more gen¬ 
eral form with three arbitrary external momenta (one 
momentum is always fixed by momentum conservation). 
In this case, also scattering processes of an electron that 
changes its patch with an electron that stays in the same 
patch contribute. As discussed in Sec. IV A, the phase 
space of these processes strongly depends on quantitative 
details of the band structure and of the patch geometry. 
From these considerations, we see that that focusing on 
the limit kpl 1 not only allows for a controlled sum¬ 
mation of a minimal set of diagrams but also renders 
the final result independent of the aforementioned mi¬ 
croscopic details. 

In the explicit summation and evaluation of the dia¬ 
grams, which can be found in the Appendix, the compo¬ 
nent of the patch momentum perpendicular to the Fermi 
surface is constrained by an energetic cutoff A. The lon¬ 
gitudinal limits of the patches is absorbed by introducing 
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the density of states pr per patch. For simplicity, we as¬ 
sume that the density of states does not depend on r 
and introduce pp = Spr, the density of states per Fermi 
surface. In Sec. V, we will also comment on the modifi¬ 
cations in case of an imbalance in the density of states 
on the different Fermi surfaces. 

Finally, we note that the two scattering matrices in 
Eqs. (34) and (35) that determine the behavior of the 
superconductor under the influence of impurity scatter¬ 
ing are made unique by the phase convention (26): Ob¬ 
viously, Srr' is already invariant under arbitrary C/(l)- 
gauge transformations of the wavefunctions, whereas Vtt' 
is fixed due to Eq. (26) since states with the same A and 
j but opposite p always come in pairs. 


V. RESULTS 

Let us begin our discussion with the clean case. Cal¬ 
culating the Ginzburg-Landau expansion in the absence 
of any disorder reproduces the results of Ref. 25: For at¬ 
tractive (repulsive) pair hopping, i.e., U < 0 (U > 0 ) in 
Eq. (30), the («■' ) superconductor is realized. In 
both cases, the transition temperature assumes the usual 
BCS form 


2e'^ 2 

T,,o = A-e (36) 

TT 

with 7 denoting Euler’s constant. 

To proceed with the impact of disorder on supercon¬ 
ductivity, we first focus on the case of nonmagnetic disor¬ 
der which is characterized by the vertex in Eq. (16). For 
attractive pair hopping, we again obtain the state. 
Within the approximations of Sec. IV C, the associated 
transition temperature is unaffected by disorder irrespec¬ 
tive of the values of the different coupling constants in the 
impurity vertex (16) . This result is a manifestation of the 
well-known Anderson theorem®^®. 

In case of repulsive pair hopping, V > 0, we still find 
that, as in the clean case, the s“' superconductor con¬ 
denses. In accordance with the usual expectation that 
unconventional superconductors are less robust against 
impurity scattering, the transition temperature is now 
reduced by disorder. Despite the complicated form of 
the vertex (16), it turns out that the impurity distri¬ 
bution only enters the Ginzburg-Landau expansion and, 
thus, the transition temperature, in form of the combined 
scattering rate 


and (35). Since the two scattered electrons either both 
switch the orbital or both keep their orbital character, 
either rgcro or T 2 CT 3 in Eq. (16) yields zero in S and V. 
Furthermore, we emphasize that the effective scattering 
strength and, hence, the transition temperature do not 
depend on the microscopic details of the impurity scat¬ 
tering, i.e., on the relative weight between different terms 
in Eq. (16). As the different coupling constants are pro¬ 
portional to the concentration of the impurities and the 
square of the matrix element of the associated local per¬ 
turbations, F in Eq. (37) measures the average overall 
strength of the disorder configuration. 

The transition temperature Tc is found from demand¬ 
ing that the corresponding eigenvalue in the disorder av¬ 
eraged Ginzburg-Landau expansion (32) changes sign. 
This leads to (see Appendix) 


(38) 

where tpo denotes the digamma function and Fc = 
e~'^Tcfij pf has been introduced for convenience. 
Eq. (38) yields the functional dependence TcjTcfi = 
/(F/Fc) of the transition temperature on the scatter¬ 
ing strength F. Mathematically, it reproduces the usual 
Abrikosov-Gorkov law® originally describing the critical 
temperature of BGS superconductors under the influence 
of spin-magnetic disorder. Analytically, one can easily 
show that /(I) = 0 , which identifies Fc as the critical 
scattering rate, as well as 

/(r/r.)^i-|;^ (39) 

for F/Fc —>■ 0. The functional behavior for general values 
of F < Fc is shown in Fig. 5(a) (together with the Ander¬ 
son theorem for s++ discussed above). In this figure, the 
scattering strength is measured in units of the reference 
critical scattering rate Fc^o- This quantity is defined by 
considering the situation where the patch wavefunctions 
would allow electrons to scatter between all patches in 
a way that does not break time-reversal symmetry (non¬ 
magnetic scattering). It corresponds to fermions with the 
same band structure but without spin and orbital degrees 
of freedom. Then Fc^o is the critical value of F = Sr.r' 
at which the s"' state is fully destroyed. To be explicit, 
one finds 


0 =— 

|V| 2 


In 




r/u 


Tc/n,, 


F 


c.O — 


2c,0 

ie^^ Pf 


(40) 


r = 7 ai+ 7"i+7 Bi+ 7B2+7i?- (37) 

Note that which describes the admixture of s-wave 
and spin-orbit disorder on the same site, does not enter 
at all. This follows from the strong orbital polarization 
of the wavefunctions in the patches (see Fig. 1), which we 
have idealized to be perfect in our calculation. It is easily 
seen upon inserting the term with 7 ^^ in Eq. (16) into 
the two relevant scattering matrices defined in Eqs. (34) 


The enhancement Fc = 4rc,o of the critical scattering 
rate of s“' as compared to the reference model is due 
to the locking between the crystal momentum and the 
spin/orbital degree of freedom. To illustrate this point, 
consider the simplest case of spin-trivial s-wave impuri¬ 
ties, i.e. only take the first term in Eq. (16). Recalling the 
spin and orbital textures of Fig. 1, an electron from, e.g., 
patch (1,2,4-) in Fig. 3 cannot scatter into one of the 
patches of the perpendicular lobes of the Fermi surfaces 
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(a) nonmagnetic scattering 
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Figure 5. The critical temperature as a function of the scattering rate F for the different superconductors in case of (a) 
nonmagnetic and (b) magnetic disorder. The scattering rates F in part (a) and (b) are defined in Eqs. (37) and (41), respectively. 
Tafi and Fc.o are the critical temperature (36) of the clean system and the critical scattering rate (40) of the reference model 
(green dashed line in (a)) defined in the main text. 


since it cannot change its orbital character. Within the 
same lobe, it can only scatter into states of the patches 
( 1 , 1 , +) and ( 2 , 2 , —) as a consequence of the spin polar¬ 
ization of the wavefunctions. We thus see that the spin 
and orbital polarization of the Fermi surfaces reduces the 
possibilities for scattering events as compared to naively 
neglecting the spin and orbital degrees of freedom of the 
electrons. Note that this mechanism of protection against 
impurity scattering can in some systems even rule out all 
scattering processes that are pair breaking as has been 
demonstrated in Ref. 16. 

By independently rescaling the intra- and interband 
scattering rate (see Appendix), we have verified that only 
the interband processes affect the transition temperature 
of the 5 “' superconductor. This is expected as, in the 
absence of interband scattering, the and s"' super¬ 
conductors have to behave identically and, hence, both 
be unaffected by disorder due to the Anderson theorem. 

Let us next discuss magnetic disorder which is, in its 
most general form, defined by the impurity vertex in 
Eq. (17). Similarly to the nonmagnetic case, the impurity 
distribution only enters in form of the single variable 

r = Tb + 7 " + iW + + TSi + IB 2 (41) 

in the expansion (32) of the average free energy. This 
shows that also the impact of magnetic scattering on 
superconductivity does not depend on microscopic de¬ 
tails of the impurities involved. It is even independent of 
whether the magnetic moments are aligned in-plane (first 
three terms in Eq. (41)) or perpendicular to it (remain¬ 
ing four coupling constants). Note that 7 ™ in Eq. (17b) 
does not contribute at all, which is, exactly as in case of 
7 ^J above, due to the spin and orbital polarization of the 
patch wavefunctions. 


We find that, as before, E < 0 (E > 0) leads to pure 
(s“' ) superconductivity. However, in case of mag¬ 
netic impurities, both superconductors are affected by 
disorder with the transition temperature again following 
from Eq. (38) where Tc = 2rc,o and Tc = dTc^o in case 
of 5 “'"'' and 5 “' , respectively. The resulting dependence 

Tc/Tcfi = /(T/Tc) of the transition temperature on the 
scattering strength T in Eq. (41) is shown in Fig. 5(b). 

Interestingly enough, the unconventional s“' super¬ 
conductor is (by a factor of 2 ) more robust against mag¬ 
netic scattering than the conventional s’*"’*' state. To un¬ 
derstand this, we have again analyzed which type of scat¬ 
tering processes act pair breaking in the two supercon¬ 
ductors. One finds that is suppressed both by inter- 
and intraband scattering, whereas the s“' state is only 
sensitive to intraband processes. This leads to the rela¬ 
tive factor of 2 mentioned above. We refer to Ref. 18 for 
another example of the stability against TRA interband 
scattering of a superconductor with an order parameter 
that changes sign between two Fermi surfaces. The ad¬ 
ditional factor of 2 with respect to Tc^o is, again, due to 
reduction of possible scattering processes resulting from 
the spin and orbital polarization of the patch wavefunc¬ 
tions. 

Finally, a few remarks are in order that refer both to 
magnetic and nonmagnetic scattering. Firstly, taking the 
spin and orbital polarization of the patch wavefunctions 
to be perfect constitutes a very good approximation as 
can be seen in Fig. 1. Introducing deviations from per¬ 
fect polarization would have the following major modifi¬ 
cations: The impact of disorder on the superconductors 
would become less universal in the sense that it would not 
be just the sum of all coupling constants as in Eqs. (37) 
and (41) determining the transition temperature. Fur- 






12 


thermore, the enhancement Fc/rc^o of the critical scat¬ 
tering rates with respect to the reference model would be 
reduced as the orbital/spin-momentum locking would be 
weaker. 

Secondly, having slightly different densities of states on 
the two Fermi surfaces would not crucially alter our con¬ 
clusions. Again, the results would become less universal, 
since, already in the clean system, we would get an ad¬ 
mixture of 5 “'"“'" and 5 “' superconductivity. Then, e.g., 
for V < 0 and nonmagnetic disorder, the order parame¬ 
ter of the resulting superconductor would have identical 
signs on both Fermi surfaces but slightly different mag¬ 
nitudes. We would find that this superconductor is de¬ 
stroyed by disorder at some finite but very large critical 
scattering strength. 

Furthermore, we emphasize that we have discussed the 
most general model of local, J-correlated nonmagnetic 
and magnetic disorder here. Any more specific model is 
thus covered by Eqs. (16) and (17). One might be inter¬ 
ested, e.g., in the case where not only the disorder distri¬ 
bution but also every single impurity has the full symme¬ 
try of the system (only 7 ) 4 ^, 7 ^ 4 ^ and 7 ™ in Eq. (16) are 
finite). To give another example, isotropic in-plane spin- 
magnetism corresponds to 7 ^ in Eq. (17a). Eqs. (37) 
and (41) show that, despite the plethora of possible mi¬ 
croscopic scenarios in the two-band model of oxide inter¬ 
faces, only the four different combinations V > 0/E < 0 
and magnetic/nonmagnetic disorder need to be discussed 
separately. 


VI. CONCLUSION 


rates, the form of the superconducting order parameter 
is unaffected by disorder. Despite the multitude of possi¬ 
ble scattering channels, the information of the impurity 
statistics only enters in form of the combined scattering 
rates F defined in Eqs. (37) and (41) for nonmagnetic 
and magnetic scattering, respectively. Our results for 
the transition temperature of the disordered supercon¬ 
ductors, which are summarized in Fig. 5, lead to the fol¬ 
lowing two central conclusions: Firstly, the critical value 
of the scattering rate F, where both superconductors in 
case of magnetic and in case of nonmagnetic impu¬ 
rities are fully destroyed, is enhanced as compared to a 
reference model that naively neglects the spin and orbital 
polarization. Although just being a factor of order 1, this 
enhancement of stability against disorder might be rele¬ 
vant for oxide interfaces where the mean-free path^®’^® is 
expected to be comparable to the coherence length^. Sec¬ 
ondly, contrary to the expectation that unconventional 
superconductivity is less robust in the presence of dis¬ 
order, we find that the critical scattering rate of s“' is 
twice as large as that of the superconductor in case 
of magnetic disorder. We emphasize that magnetic im¬ 
purities are expected to be particularly relevant in oxide 
interfaces as the system has the tendency towards the 
formation of a magnetic instability®'^^®®. For example, 
oxygen vacancies have been shown^^^i ^ induce local 
magnetic moments on the Ti sites. Taken together, these 
findings indicate that the presence of superconductivity 
despite the coherence length and mean-free path being 
comparable does not allow for ruling out unconventional 
pairing. 


To summarize, we have considered a model for the itin¬ 
erant electrons at oxide interfaces that includes the Ti 
Sdxz and 3dyz orbitals which are expected to be most 
relevant for describing superconductivity ^®’^°’®®“®®. De¬ 
pending on microscopic details of the electron-electron 
interaction, this model yields^® a conventional or an 
unconventional s® superconductor. 

To describe disorder scattering, we have taken the most 
general ^-correlated impurity distribution that is consis¬ 
tent with the spatial symmetries of the system (point 
group Ciy and 2D translation) as well as with time- 
reversal symmetry. Note that these symmetries only have 
to be restored on average such that a single impurity can 
locally break any symmetry of the system. The combina¬ 
tion of spin and orbital degrees of freedom leads to many 
different possible scattering channels as can be seen in the 
impurity vertices in Eqs. (16) and (17) for nonmagnetic 
(TRS) and magnetic (TRA) impurities, respectively. We 
have discussed how these different processes can emerge 
due to vacancies or the formation of local magnetic mo¬ 
ments. 

As very different behavior is expected, we have stud¬ 
ied the impact of magnetic and nonmagnetic disorder 
on the two superconductors s"*'"'" and s"' separately. 
In all cases, we find that, below the critical scattering 


Appendix: Summation of the diagrams of the 
Ginzburg-Landau expansion 


This appendix provides details on the calculation of 
the coefficients dxy{T) in the disorder-averaged expan¬ 
sion of the free energy in terms of the order parameters as 
defined in Eqs. (32) and (33), and diagrammatically rep¬ 
resented in Fig. 4. The approach we present is tailored 
to the description of oxide heterostructures in the patch 
approximation (Sec. IV A) of the Fermi surfaces. How¬ 
ever, the described procedure can be modified straight¬ 
forwardly for application to other multiband systems as 
long as the Fermi surfaces can be partitioned in a similar 
way that allows for a formulation of the problem in terms 
of scattering rates that are constant within each patch. 

In the presence of disorder, the full propagator (cf. 
Fig. 4(b)) in patch r acquires a finite lifetime. It is mod¬ 
ified by a finite self energy Et-(p, a;„). 


Grip,UJn) 


1 


(A.l) 


where €r{p) denotes the disperson in patch r. We con¬ 
sider weak disorder such that diagrams with crossed im¬ 
purity lines are suppressed by a factor {kpl)~^ and calcu¬ 
lating the self energy in Born approximation is sufficient. 
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These assumptions result in the expression 


Sr(P,W„) = -iTTS\g-D.{uJn)'^PT'SrT' (A.2) 


for the self energy in patch r, where p-r is the density of 
states in patch r, and the corresponding scattering ma¬ 
trix Srr' of scattering processes contributing to the self 
energy in the patch approximation is defined in Eq. (34). 
As introduced in Sec. IV C, we assume that the density 
of states does not depend on r and use pp = in the 
remainder of this appendix. 


The expansion coefficients d\\i (T) in the presence of 
disorder contain contributions of all possible noncross¬ 
ing combinations of the allowed scattering processes of 


Kramers partners (cf. the definition of Vtt' in Eq. (35)), 



(A.3) 

For within the patch approximation, the impurity ver¬ 
tex is a constant Vrr' for all scattering processes of 
Kramers partners from patch r to patch r', the momen¬ 
tum integrals in Eq. (A.3) decouple and it is useful to 
define 

^ ^ (P; ^n)f^TK ( Pi ^n) (A.4) 

P 

which evaluates to 

/ X __ '^Pf _ 

" 8|a;„| + ^ sign(u;„) ’ 

(A.5) 

where we used er^i—p) = ^t[p)- 

Exploiting the patch approximation and following the 
notation introduced in Sec. IV C, the coefficients in the 
Ginzburg-Landau expansion can be written in the form 


d\\' 





T 

2 


EEE dV 0(A,j,7)) (Wn)V(A J,r/)(A'j'.r;')C(A'(Wn) 

j,j' V,V' 


T 

2 


EEE E 

“n A" j,j' ,j" ri,ri' ,ri" 


(A.6) 


r 


suggesting to introduce, analogously to S and V, nested 
diagonal matrices C containing the propagator products 


Ctt' — drr' ■ 


(A.7) 


With these definitions at hand, we can conveniently ex¬ 
press Eq. (A. 6 ) in matrix form. 


C^AA' =-|EEEw' 


j,j' 71,7]' 


^ (CV)™C 


m—0 


J {X,j,ri){X',j’,7]') 

(A.8) 

and therefore, the results for the coefficients d\\i(T) can 
be obtained from matrix multiplication and matrix in¬ 
version. 


= (i-cvr'c 


AaA' 


jj' v,v' 


- (A,i,J))(A'j', 77 ') 

(A.9) 


This result can be staightforwardly applied to other 
multiband superconductors, if an analogous partitioning 
of the Fermi surface allows to formulate the problem in 
terms of scattering rates that are only characterized by 
respective patch indices and, apart from that, do not de¬ 
pend on the momenta of the scattered states. 

We investigate the different roles of intraband and in¬ 
terband scattering processes for nonmagnetic and mag¬ 
netic disorder by artificially discriminating intraband 
(A = A') and interband (A 7 ^ A') scattering processes 
by setting 


5, 


(Aj'.»))(A'j',?)') 


'^{X,j.rj){X',j',r}') 


{ Xintra‘S(X,j,ri)(X' ,j' ,ri') 1 A — A , 
XinteT<S(XJ,?j)(X' J' ,ri') 1 A yf A , 

(A.IO) 

{ XintraV(A,j,r/)(A',j',7)')i A = A , 
XinterV(Aj',r;)(A'j'.r;') ) A A . 

(A.Il) 
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In order to reproduce Eq. (38) and the dependence of 
the transition temperature Tc on the scattering rate F as 
shown in Fig. 5, the results discussed in the remainder 
of this appendix have of course to be understood in the 
sense that Xintra “ Xinter — 1- 

The resulting diagrams in the patch approximation 
satisfy dii{T) = d 22 {T) and di 2 (r) = d 2 i{T), as well as 
dll € M and di 2 G K, and thereby, the free energy (32) 
can be diagonalized as 


r-l 


(•^)dis + di 2 iT)\ ^ (Ai + A 2 ) 


— + dii(T) - di2(r) ^ (Ai - A 2 ) 


(A.12) 


presence of nonmagnetic disorder is 


0 = 


1^1 


Pf 

2 L 




1 I PFXinterF^ 




(A. 17) 


meaning that only interband scattering processes affect 
the transition temperature of the s“' pairing state, 
and the suppression of the transition temperature Tc 
with increasing scattering rate F is described by an 
Abrikosov-Gorkov law® with critical scattering rate F^ = 
e-~^Tcfi/P f- 

h. Magnetic disorder In the presence of time- 
reversal antisymmetric (TRA), i.e. magnetic, disorder, 
we find for the expansion coefficients 


The resulting transition temperature for either s’'""'" 
(Ai = A 2 , in case of V < 0) or s’'"” (Ai = —A 2 , in case 
of V >0) superconductivity is found from demanding 
that the corresponding eigenvalue change sign at T = Tc- 
a. Nonmagnetic disorder In case of disorder poten¬ 
tials that are symmetric under time-reversal (TRS), that 
is, for nonmagnetic disorder, we find 


dii(T) 

di2(T) 




di^XinterF \ 

TT /. 


(A.13) 


(A.14) 


and consequently, the transition temperature for the s’'”'’ 
superconducting state in the presence of nonmagnetic im¬ 
purities is given by the condition 




(A.15) 


resulting in a constant transition temperature that coin¬ 
cides with the transition temperature in the clean case. 


2e'> 2 

Tc = A-e ^ = Tc,o, 

TT 


(A.16) 


rfii(T) 

di2{T) 


Pf 

' 4 L' 
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4T ) 


V'o ( 
V'o ( 
V’o ( 


Pf ' 
' 4 . 


1 PF(Xintra + Xinter)F'^ 

2 

1 
2 


4T 

/^i^XintraT \ 

TT J 


1 PF(Xintra + Xinter)FN 

2 4T J 


(A. 18) 


(A. 19) 


The transition temperature of the s^"'" pairing state is 
thus determined by the condition 


0 =— 

|R| 2 




1 PF(Xintra T Xinter)FN 


4T 


(A.20) 

yielding an Abrikosov-Gorkov-type suppression with a 
critical scattering rate of Fc = e~^Tcfl/{2pF) since 
s"'"^ superconductivity in the presence of magnetic dis¬ 
order is destroyed by intraband as well as by interband 
scattering processes. In contrast, s"' superconductivity 
is only affected by intraband scattering processes, as seen 
from the corresponding condition 


0 =— 

|R| 2 L 




P^XintraT N 


(A.2I) 


that is, we reproduce the Anderson theorem® ®. The cor- and thus the respective critical scattering rate is again 
responding condition for s"' superconductivity in the given by Fc = e~'^Tcfijpp- 
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